Protein crystallization in confined geometries 
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We studied the crystallization of a globular protein, lysozyme, in the cubic phase of the lipid 
monoolein. The solubility of lysozyme in salt solution decreased by a factor of ~ 4 when confined 
in cubic phase. Calculations and Monte Carlo simulations show that this can be explained by the 
confinement of lysozyme molecules to the narrow water cells in the cubic phase. 



We report a study of the crystallization of lysozyme, 
a globular protein, in the cubic phase of monoolein, an 
uncharged lipid. The protein solubility was found to de- 
crease by a factor of ^ 4. We account for this obser- 
vation by the confinement of lysozyme molecules within 
nanoscopic water 'cells' in the cubic phase. The result- 
ing entropy-driven rise in their chemical potential can ex- 
plain our observation of solubility suppression compared 
to the lipid-free case. The generic nature of this expla- 
nation means that our results are of potential interest 
to three rather different audiences. First, and most gen- 
erally, scientists interested in the basic physics of con- 
finement P] may find the lysozyme-monoolein system a 
useful model for detailed study. Secondly, there is grow- 
ing interest in directing the self-assembly of synthetic or 
natural nanoparticles using various means of confinement 
(e.g. at an interface |3| or in emulsion droplets Our 
work suggests that lipid cubic phases can be useful in 
this regard. Thirdly, and most specifically, lipid cubic 
phases offer an environment for globular protein crystal- 
lization Q. (Their use for crystallizing membrane pro- 
teins 1^, which reside in the lipid bilayers, relies on dif- 
ferent physics 

Lysozyme crystallization in lipid cubic phase has been 
reported before 0, although to date there is no com- 
plete phase diagram or physical explanation of the mech- 
anism. We have studied the phase behavior of lysozyme 
solution as a function of protein and salt concentration 
in a monoolein cubic phase, and compared it with the 
case without lipids. The cubic phase was found to lower 
the solubility dramatically. This result was then explored 
using simple calculations and Monte Carlo simulations. 

Lysozyme and monoolein (MO) were used as purchased 
from Sigma- Aldrich. Protein concentration in 0.05 M 
sodium acetate buffer was determined, after filtration 
through a 0.1 /xm Millipore filter, by UV spectropho- 
tometry (specific absorbance v4280nm = 2.64 ml/mg.cm). 
NaCl solutions of known concentration were then added. 

We first determined the bulk solubility boundary of 
lysozyme/salt without lipid from the dissolution of 'seed' 
crystals which occurs at lower protein and salt con- 
centrations than the bulk nucleation boundary (where 
crystals are first seen) due to the nucleation barrier 
Next protein/salt solutions with concentrations be- 
low the bulk solubility boundary was added to and vig- 



FIG. 1: An optical micrograph of lysozyme crystals in 
monoolein cubic phase. Note the curved piece of solid in- 
clusion; other crystals grow next to surfaces. 



orously mixed with MO powder melted on a microscope 
slide. A highly viscous cubic phase was observed to form 
almost instantly. The excess solution was then removed 
and the cubic phase sealed under a cover slip (sample 
thickness 100 /xm) for observation in an optical micro- 
scope at 18°C. All samples were transparent and optically 
isotropic under cross polarizers. Within two months, 
crystals were observed in many samples, e.g. Fig. ^ 

Fig. [5] summarizes our observations in the form of a 
phase diagram. We give lysozyme concentrations, c, in 
mg/ml, and in terms of an effective volume fraction, (j), 
by modelling the globular protein as a sphere of diameter 
3.4 nm (which reproduces the molecular volume). Note 
that for the data points recording behavior in cubo, the 
protein/salt concentrations given are those of the parent 
solution into which MO powder was added. The bound- 
ary above which we observed crystals in cubo (a kinetic 
nucleation boundary) lies entirely below the bulk solu- 
bility boundary. Since we expect the in cubo solubility 
boundary to lie at even lower protein/salt concentrations, 
our main experimental finding is that lipid cubic phase 
significantly depresses lysozyme solubility. 

To understand this observation, we first note that the 
bicontinuous water in the cubic phase is made up of cells 
connected by narrow 'necks'. The presence of lysozyme 
changes the MO cubic phase structure from PnSm to 
Im3m |lO|| . a schematic projection of which is shown in 
Fig. inja) with an inter-cell distance of 12.5 nm ^3 ^tnd 
bilayer thickness of 3.5 nm We have also shown a 
number of lysozyme molecules (hatched circles) modelled 
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FIG. 2: The phase diagram of lysozyme/NaCl in bulk solu- 
tion and in MO cubic phase. Filled triangles (a) represent the 
solubility limit in bulk solution. Other symbols represent be- 
havior in the cubic phase: homogeneous solution (x), crystals 
(□), and other precipitates (A). The lines are guides to the 
eye. Inset (same axes), the concentration of the samples in 
the cubic phase is multiplied by 4, to nearly superimpose the 
nucleation boundary in cubo on the bulk solubility boundary. 



as spheres with diameter cr = 3.4 nm. Since the bending 
rigidity of the bilayers is fee ~ ^ksT at room temperature 
proteins cannot move between cells without incur- 
ring significant cost in distorting the connecting 'necks'. 
Thus the lysozyme will largely be confined in nanoscopic 
water 'cells' with linear dimension < 10 nm , ea ch of which 
can hold up to 7 lysozyme molecules [l^ . The cen- 
ter of a lysozyme molecule does not come closer than 
^ (t/2 = 1.7 nm to a bilayer. While such excluded vol- 
ume is a small effect in the bulk, it is highly significant 
inside a water 'cell' in the cubic phase, leading to a large 
increase in the effective protein concentration. This ex- 
plains why crystallization occurred at lower average pro- 
tein concentration. 

Enquiring where the crystals actually nucleate gives 
further understanding. The critical nucleus in our exper- 
iments is expected to consist of > 10 lysozymes 0|, more 
than an undistorted water 'cell' can hold. Moreover, a 
(tetragonal) unit cell contains 8 molecules and it is 
unlikely that crystals can nucleate in pores smaller than 
a few unit cells Experimentally, we often observed 
crystals next to solid inclusions or next to surfaces, Fig.^ 
We therefore propose a role for defects. 

The picture is as follows. The protein in bulk so- 
lution has a concentration-dependent chemical poten- 
tial /i(c). Confinement in cubo increases the effective 
lysozyme concentration to c' > c, and its chemical po- 
tential to /i(c') > /i(c). There are defects in the cubic 
phase of mesoscopic sizes (^ lysozyme diameters), due 



to grain boundaries, distortion by inclusions, or surfaces. 
In these defects, the chemical potential of lysozyme is ini- 
tially at /x(c), but diffusion from the 'loaded' cubic phase 
rapidly increases this to close to /i(c'). Crystals nucleate 
from this concentrated solution residing in the defects. 
Subsequently, as a crystal grows, it will break the sur- 
rounding matrix; the elasticity of the cubic phase may 
therefore affect the rate of crystal growth. 

Support for the above picture comes from mea- 
suring the time-dependent protein concentration in 
lysozyme/salt solution in contact with the cubic phase 
formed by mixing it with MO. We observed (data not 
shown) a rapid increase in concentration that continued 
for many hours, showing that indeed //(c') > //(c). 

To explore quantitatively our proposed picture, we 
used Monte Carlo (MC) simulations to study attractive 
hard spheres confined within a 'model cubic phase'. To 
do this, we first need a numerical estimate of the excluded 
volume effect. The center of a lysozyme molecule con- 
fined to a water 'cell' in the cubic phase cannot explore 
the whole of the cell volume, Vc, but only the accessible 
volume, Va, given by taking off the volume of a layer of 
thickness a/2 from Vc- The ratio a = Vc/Va characterizes 
the excluded- volume effect. 

Estimating a for the ImSm structure involves non- 
trivial geometry and knowing the exact lattice parameter, 
a, pertaining to our experiments. The value of a used in 
Fig. 13 is for fully- hydrated MO cubic phase equilibrated 
at 25*^0 These conditions almost certainly do not 

apply (after mixing at 25° C we immediately removed ex- 
cess solution and observed at 18°C). Nevertheless, using 
a — 12.5 nm [loj as an estimate, we can approximate 
the accessible volume in each water 'cell' as a regular 
square bipyramid of side 8.8 nm, Fig.|2Ib). The centres 
of lysozyme molecules can explore a smaller bipyramid of 
side 8.8 - cr = 5.4 nm, giving aimSm ~ (8.8/5.4)^ = 4.3 
[l^ . Literature data for MO under different conditions 
without ^3 and with lysozyme suggest that we can 
expect variations of up to ^ ±30% in the lattice param- 
eter, so that a ^ 3 — 11. 

In our simulations, we used small spheres to model 
lysozyme, while 125 larger spheres fixed on a primitive 
cubic lattice modelled the geometric confinement effect of 
the MO cubic phase. Such a model has been used before 
to study fluids in porous media 0|. The small sphere 
diameter and large-sphere simple cubic lattice constant 
were taken to be 3.4 nm and 12.5 nm respectively, to 
agree with the known lysozyme molecular volume and 
the cubic phase with lysozyme lattice constant at 25'^C 
[lo| . We simulated a model with a particular a by using 
a value for the large sphere diameter such that the small- 
sphere centers can access a^^ of the total space between 
large spheres. Periodic boundary conditions apply. 

The small spheres (diameter a) interact via a hard core 



FIG. 3: (a) A 2-d projection of the ImSm cubic phase, with 
'necks' (a-d) connecting water 'cells'. Annular regions are 
lipid bilayers (with lipids shown schematically). Shaded cir- 
cles are cross sections of 'necks' in and out of the plane of the 
paper. Hatched particles are (arbitrarily-placed) lysozymes 
with diameter a — 3.4 nm. The lattice parameter is 12.5 nm. 
The larger square (broken line) is the base of a regular square 
bipyramid used to model a water 'cell'. The smaller square is 
the base of the (smaller) bipyramid accessible to the centres 
of lysozyme molecules, (b) A bipyramid 'cell' with lysozyme 
molecules. 



and a square-well attraction: 



-oo, r < a 
u{r) = { e a <r < {1 + X)a 
r>(l + A)(T 



(1) 



where A and e are the range and depth of the square well 
in units of a and fc^T respectively. We took A = 0.08 [l9|. 
Increasing the salt concentration experimentally corre- 
sponds to increasing e; in our range of salinity, we expect 
< e < 3 10 . The small spheres behave as hard spheres 
to the fixed large spheres. 

We investigated the coexistence of a fluid of such at- 
tractive small hard spheres confined by a fixed array of 
large hard spheres with an unconfined crystal of the at- 
tractive hard spheres. This is the thermodynamic situ- 
ation of crystals nucleating in mesoscopic defects filled 
with concentrated protein solution 'fed' from the cubic 
phase. Rather than simulating the free energies, we 
used an approximation to calculate coexistence by simply 
equating chemical potentials and neglecting the equal- 
pressure condition (see for a careful justification). 

First, in the spirit of the cell model (ij we can approx- 
imate the volume accessible to the centre of a particle in 
the crystal to be nX^a^ /6, so that the chemical potential 
of the crystal can be written as, 



Ms = Mo - ns 



— — SknT In X, 
2/ 



(2) 



where Us is the number of particles whose centers lie 
in the region a < r < (l-t-A)cr from a given parti- 
cle. Here, /io — ksT \n{6A^ /ira^) where A is the de 
Broglie (thermal) wavelength. Note that with this choice 
of fiQ, the ideal-solution chemical potential is given by 
Mi = Mo + ksTliKl). 
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FIG. 4: The excess chemical potential for a fluid with A = 0.08 
at e = 1.5: a = 1 (unconfined, □), 2 (A), 4 (o), 6 (V), and 8 
(o) from simulation. The lines are guides to the eye. 



Next, consider the fluid of confined small spheres. The 
leading order correction to the ideal limit is 



fif = kBTlii{paA^) = Ml + ksTlna , 



(3) 



where pa is the number density of small spheres in the vol- 
ume accessible to them, and a^^ is the ratio of this acces- 
sible volume to the total volume between large spheres. 
In our simulations, a = 2.0 j4.0 , 6.0 or 8.0. We used the 
'Widom insertion method' [23| to obtain Afj, = — fj,i 
at different values of e with and without confinement, 
Fig. 21 Note that the confined results indeed go to 
Ap/ksT = In a in all cases when (/) — > 0. 

Now we calculate fluid-crystal coexistence at each 
value of e by equating the crystal chemical potential. 
Equation 121 with Ug — 12, to the simulated fluid chem- 
ical potential j20|]. The low-density branches with and 
without conflnement are shown in Fig. [S| Increasingly 
pronounced confinement indeed systematically decreases 
the solubility in a way very reminiscent of the data shown 
in Fig. 13 

Quantitatively, we show in the inset that multiplying 
the volume fraction ^ by a leads to overlap of the phase 
boundaries for a(f> < 0.2. Note that the same result is 
obtained (data not shown) if we map e to the salt con- 
centration via measured second virial coefficients [l^l23| . 

This result suggests that we should attempt the same 
procedure on the experimental data. This is problem- 
atic. The solubility boundary in cubo cannot be mea- 
sured: Fig. 121 shows a nucleation boundary. If we nev- 
ertheless proceed, we find that multiplying the lysozyme 
concentration by a factor of 4 shifts the confined nucle- 
ation boundary to the bulk solubility boundary (inset. 
Fig- El- The in cubo solubility boundary will undoubt- 
edly occur at lower concentrations than the nucleation 
boundary. We therefore expect that a 'shift factor' of 
> 4 will be needed to translate the (unavailable) in cubo 
solubility boundary to overlap with the bulk solubility 
boundary |24l| . Given that we expect a ~ 3 — 11 in our 
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FIG. 5: The simulated equilibrium crystallization boundary 
with a = 1 (unconfined, □), 2 (A), 4 (o), 6 (V), and 8 (o). 
The lines are guides to the eye. In the inset, the volume 
fraction of each system is multiplied by a. This superimposes 
all five boundaries at ad) < 0.2. 



cubic phase, this expectation is not inconsistent with our 
simulation results. 

In reporting our experimental findings, we have so far 
not mentioned any time-dependent features. Two aspects 
are worthy of note. First, we often observed amorphous 
precipitates before the appearance of crystals. Unlike 
crystals, these precipitates would occur rather uniformly 
throughout a sample, and often appeared as 'stars' in 
the microscope. Sometimes, growing crystals would 'con- 
sume' surrounding amorphous precipitates. The origin 
and nature of these precipitates are unclear. Secondly, 
although it is difficult to be quantitative, it is our im- 
pression that the in cubo growth of crystals is a slower 
process compared to that in bulk solution. This may, of 
course, be a reflection of the fact that 'feeding' lysozyme 
into the growth medium residing in defects is a slow pro- 
cess. 

To conclude, we have studied the phase behavior of 
lysozyme/salt in MO cubic phase. The latter lowers the 
solubility of the protein. This is predominantly a confine- 
ment effect: experiments and simulations suggest that a 
single 'excluded volume scaling factor' can shift the solu- 
bility boundary back to its bulk position. The altered ki- 
netics of crystallization in cubo, relying as we suggest on 
defects and other properties of the cubic phase, may of- 
fer a new means of controlling protein crystal nucleation 
and growth. The generic, entropic, nature of our pro- 
posed explanation suggests that lipid cubic phases may 
be useful for manipulating a range of nanoparticles, and 
for studying the basic statistical physics of confinement. 
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